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Abstract. Local energy error estimates for the finite element method for el- 
liptic problems were originally proved in 1974 by Nitsche and Schatz. These 
estimates show that the local energy error may be bounded by a local ap- 
proximation term, plus a global "pollution" term that measures the influence 
of solution quality from outside the domain of interest and is heuristically of 
higher order. However, the original analysis of Nitsche and Schatz is restricted 
to quasi-uniform grids. We present local a priori energy estimates that are 
valid on shape regular grids, an assumption which allows for highly graded 
meshes and which much more closely matches the typical practical situation. 
Our chief technical innovation is an improved superapproximation result. 



1. Introduction 

In this note we prove local energy error estimates for the finite element method 
for second-order linear elliptic problems on highly refined triangulations. Most a 
priori error analyses for the finite element method in norms other than the global 
energy norm place severe restrictions on the mesh. In particular, such error analyses 
are most often carried out under the assumption that the grid is quasi uniform, that 
is, all simplices in the mesh are required to have diameter equivalent to some fixed 
parameter h. The typical practical situation is rather different. Many (especially 
adaptive) finite element codes enforce only shape regularity of elements, meaning 
that all elements in the mesh must have bounded aspect ratio. Though it places a 
weak restriction upon the rate with which the diameters of elements in the mesh 
may change, shape regularity allows for the locally refined meshes that are needed 
to resolve the singularities and other sharp local variations of the solution that 
occur in the majority of practical applications. 

In the work |NS74| of Nitsche and Schatz, local energy error estimates were 
established for interior subdomains under the assumption that the finite element 
grid is quasi-uniform. Such local energy estimates are helpful in understanding 
basic error behavior, especially "pollution effects" of global solution properties on 
local approximation quality, and they also provide an important technical tool in 
many proofs of pointwise bounds for the finite element method (cf. [SW95J). In 
addition, the most relevant error notion in applications is often related to some 
local norm or functional instead of to the global energy error, as evidenced by the 
recent surge of interest in ensuring control of the error in calculating "quantities 
of interest" in adaptive finite element calculations instead of merely controlling the 
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default global energy error (cf. [BR01 ). As a final example of the applicability 
of local energy estimates, we mention that the estimates of |NS74j have been used 
to justify certain approaches to parallelization and adaptive meshing (cf. [BHOOJ). 
Thus local energy estimates are of broad and fundamental importance in finite 
element theory. 

Here we prove local energy error estimates under the assumption that the fi- 
nite element triangulation is shape regular instead of under the more restrictive 
assumption of quasi uniformity required in |NS74j . In other words, we essentially 
prove that the results of Nitsche and Schatz hold under the restrictions typically 
placed upon meshes in practical codes, which in particular allow for highly graded 
grids. Our main innovation is a novel "superapproximation" result which we state 
and prove in §2. In §3 we then prove a local energy bound that is valid on grids 
that are only assumed to be shape- regular. As in [NS74j . our results are valid for 
operators that are only locally elliptic, so that the PDE under consideration may be 
degenerate or change type outside of the domain of interest. In contrast to [NS74J, 
the results we present here are valid up to the domain boundary, allow for nonho- 
mogeneous Neumann, Dirichlet, and mixed boundary conditions, and also require 
only Lqo regularity of the coefficients of the differential operator. 

2. An improved superapproximation result 

An essential feature of the proofs of local error estimates given in |NS74] , and 
also of essentially all published proofs of local and maximum-norm a priori error 
estimates for finite element methods, is the use of superapproximation properties. 
In essence, superapproximation bounds establish that a function in the finite ele- 
ment space multiplied by any smooth function can be approximated exceptionally 
well by the finite element space. 

In order to fix thoughts, we shall in this section assume for simplicity that f2 C M. n 
is a polyhedral domain; a more general situation is considered in §3 below. Let Th 
be a simplicial decomposition of f2. Denote by Ht the diameter of the element 
T G Th. We assume throughout that the elements in Th are shape-regular, that is, 
each simplex T G Th contains a ball of diameter c\ hx and is contained in a ball 
of radius C\hr, where c\ and C\ are fixed. Let also S r h be a standard Lagrange 
finite element space consisting of continuous piecewise polynomials of degree r — 1. 
We shall use standard notation for Sobolev spaces, norms, and seminorms, e.g., 
\\Am m = (Uu 2 + |V U | 2 ) Axf/\ \u\ wm = (£ w=fc ||^||^ (0) )^, etc. 

A standard superapproximation result is as follows. Let to G C°°(fi) with 
\u\ w o ( ) < Cd~\ < j < r. Then for each x £ S^, there exists rj G such 
that for each T G Th satisfying d > hx, 

(2.1) \\u X - r,\\ H t (T) < C(^\\V X \\l 2 (t) + ^llxlk(T)). 
Our modified result follows (cf. Guz06j). 

Theorem 2.1. Let to G C°°(ft) with \^\ wL{n) < Cd' j for < j < r. Then for 
each ^ £ there exists rj G SI such that for each T G Th satisfying d>hx, 

(2.2) ||c 2 x - r,\\ HHT) < C(^||V( WX )||l 2 (t) + ^llxlk(T)). 

Remark 2.2. There are two differences between (|2.1[) and (|2.2[) . First, in (|2.ip we 
consider approximation of u>x, whereas in (|2.2j) we consider approximation of oj 2 X - 



LOCAL ENERGY ESTIMATES ON SHARPLY VARYING GRIDS 



3 



Secondly, in (J2JJ) the norms on the right hand side involve only x, whereas in ()2.2j) 
the H 1 seminorm involves u>x- If we think of uj as a cutoff function, this distinction 
becomes vitally important: cox nas the same support as w 2 x, whereas the support 
of x is generally larger than that of u>x- This seemingly minor difference will allow 
us to establish local energy estimates on grids that are only assumed to be shape 
regular. 

Proof. Let Ih : C°(£l) — > be the standard Lagrange interpolant. We shall choose 
i] = Ih{(jJ 2 x) m (|2.2p . For T S 2J,, we may use standard approximation theory (cf. 
|BS02j ) to calculate 

ll^ 2 X - h(u 2 x)\\m(T) <Ch™ /2 \\uj 2 x - Ih(^ 2 X)\\w^(T) 

' ^ru n / 2 + r - 1 \ 2 i 

<Ch T < \uj X\W^(T)- 

Noting that D a x = for all multiindices a with |a| = r, recalling that ^j- < 1, 
and employing inverse estimates, we compute 

r 

C^ /2+r -Vxl«- (X) < C(^^~V 2 |^(t))IIxIIl 2( t) 

(2 .4) +C^ 2+ ^ 1 2 \\D^D^x\\lut) 

|a|=l,|,3|=r-l 

<c^||xIU 2( t) + c^" 1 £ ||^ 2 ^ X || WT) . 

|a| = l,|/3|=r-l 

We next consider the terms ||D"u; 2 D' 3 x||i oo (r) above. Since \a\ = 1, we have 
D a uj 2 = 2ujD a Lo. Let lu = f T udx so that \\uj — wHl^t) < C^tMw^t) < 
C^j-. Employing inverse estimates, we thus have 

\a\ = l,\(3\=r-l 

|/3|=r-l 

(2.5) ^Crt?/ 2 ^ 1 £ (IK«-w)P /J xlU 00 (T) + ||^ /J xlU.(T)) 

|/3|=r-l 

<C(^jIIxIIl 2 (T) + ^Pxlffi(T)) 

<C(-^2-||xIIl 2 (t) + -t\{& - w)x|ffi( T ) + — |wx| H i (T) ). 
Using an inverse inequality, we find that 

— - uo)x\hht) <-j(Hw^(t)\\x\\l 2 (t) + \\u-u\\ Lx{T )\x\hi(t)) 

(2-6) <c^(i|| x || £a(T) + ^l| x |^ (T) ) 

<C^\\x\Wt). 

Inserting (|2.6[) into (|2.5[) and the result into (|2.4[) and (|2.3p completes the proof of 
(I2~2l). □ 
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3. Local H 1 estimates 

In this section we state and prove a local H 1 estimate that is valid on highly 
graded grids. We now let be a domain in R n , and let f^o be a bounded subdomain 
of f2. We decompose dilnd^lo (if it is nonempty) into a Dirichlet portion Yd and a 
Neumann portion Tn- For the sake of simplicity, we assume that is polyhedral 
and that Tn is either polyhedral or Lipschitz. Let u satisfy 

— div(AVu) + b ■ Vu + cu =/ in Oo, 
^3 js u =g D on T D , 

du 

- — =g N on T N . 

on A 

Here A is an n x n coefhcient matrix that is uniformly bounded and positive definite 
in f2, b G ioo(r2o)™, c G Loalflo), and is the conormal derivative with respect 
to A. We also assume that SI C K™. Note that we make no assumptions about the 
differential equation solved by u outside of Q . 

Let H D0 {Q ) — {u G fl" x (fi ) : "|r D = 0}, and let H D (n Q ) = u G ff^fto) = 
u \r D = Also H<{B) — i u e -ff 1 (^o) : u|n\s = 0} for subsets B of fio- 

Thus functions in H\{B) are zero on dB \ <9fi, but may be nonzero on portions of 
dB coinciding with dft, or put in other terms, functions in H^(B) are compactly 
supported in B modulo dfl. Rewriting (|3.ip in its weak form, we find that u G 
H D (Qo) satisfies 



(3.2) 



L{u, v) :— / (AVu\7v + b ■ Wuv + cuv) dx 
Jn 

= [ fvdx- [ g N vda, v € H Dt0 (fi)nff x (O ) 
Jn Jtm 



Following |NS74] . we do not assume that L is coercive over H 1 (Q ), but rather 
we make a local coercivity assumption: 

Rl: Local coercivity. There exists a constant do > such that if B is the intersection 
of any open sphere of diameter d < do with Jlo, then L is coercive over H<(B), 
that is, for some constant C\ > 0, 

(3.3) (Cxr^uf^ < L{u,u) < CiN& 1(B)) u G H^B). 

Remark 3.1. Rl may be satisfied in one of two ways. It may happen that L is 
coercive over H (Oo) 5 i n which case no further argument is needed. Rl so long as 
a Poincare inequality 

(3-4) Hk( B) < Cd\\u\\ m{B) 

holds for balls B as in Rl having small enough diameter (cf. Remark 1.2 of |NS74| ). 
Such Poincare inequalities always hold for interior balls. If B is the nontrivial 
intersection of an open ball with Q, then (|3.4[) holds for d < d\ small enough under 
the restrictions we have placed on dfl (~l 9f2o; here d\ depends on the properties of 

dnndn . 

Next we make assumptions concerning the finite element approximation Uh of u. 
Let 7o be a triangulation such that fio C Dt£T T and T n fio $ f° r au T G T a . 
Let hx — diam(T) for T G 7~o- We denote our trial finite element space by Sd- 
We do not assume that Sd C H d (Q). In addition, we let Sd,o = Sd H -ffjj (^o) 
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be our trial finite element space. We assume that ut is the local finite element 
approximation to u on fio, that is, uh G Sd and 

(3.5) L(u - u h ,v h ) = for all v h G S D .o D fl^(fio). 

We do not explicitly fix Uh on the Dirichlet portion of the boundary, but rather im- 
plicitly assume that Uh\v D is set equal to some appropriate interpolant or projection 
of g D - 

Next we state properties that Sd and Sd,o must possess in order to prove the 
desired local energy error estimate. Let d < do be a fixed parameter, and let G\ 
and G be arbitrary subsets of Oo with G\ C G and dist(Gi,dG \ <90) = d > 0. 
Then the following are assumed to hold: 

Al: Local interpolant. There exists a local interpolant / such that for each u G 
#<(Gi), /ueS D n #<(G), and for each u G F^ i0 (fi Q ), iu G S A0 - 

^4^: Inverse properties. For each xS'S'd^G'^Ml^P^'i'i; oo, and < ^ < 
s <r with r sufficiently small, 

(3.6) llxlk|(T) < c4" S+? " ? |jx||^(T)- 

A3: Super approximation. Let w 6 G°°(O ) ni/<(Gi) with , n s < Crf - - 7 for 

integers < j < r with r sufficiently large. For each \ € <Sd,o and for each T £7~h 
satisfying d < hr, 

(3.7) \\uj 2 x - I(oj 2 x )\\ m(T) < C(^[[V(wx)|| £a (T) + ^IMU 2 (t)), 
where the interpolant 7 is as in Al above. 

Remark 3.2. Al, A2, and A3 are satisfied by standard finite element spaces defined 
on shape-regular triangular grids. Al also essentially requires that the finite element 
mesh resolve G \ Gi, i.e., that d > if maxmg^j hr with K large enough. 

We begin by proving a Caccioppoli-type estimate for "discrete harmonic" func- 
tions. Such a statement was also proved in |NS74| as a preliminary to local energy 
estimates, though the proof we give below more closely follows |SW77j . 

Lemma 3.3. Let Go C G C Flo be given, and let dist(Go,dG \ dfl) — d with 
d < 2do where do is the parameter defined in the assumption Rl. Let also Al, A2, 
and A3 hold with d = 4, and assume that uu G Sd,o satisfies 

(3.8) L(u h ,v h ) = for all v h G S D ,o D #<(fi )- 
In addition let vciS2tTc\G=jt<li -p < \. Then 

(3-9) IK|| H i(Go) <C^\\u h \\ L2{G) . 

Here C depends only on the constants in i3. 6\) and \3. 7| ) and the coefficients of L. 

Proof. We assume that Go is the intersection of a ball Bd of radius 4 with O ; the 
general case may be proved using a covering argument. Let then Gi and Gi be the 
intersections with Oo of balls having the same center as Go and having radii 4 and 
respectively, and without loss of generality let G be the corresponding ball of 
radius d. Let then ui G C Q y3 (Gi) be a cutoff function which is 1 on Go and which 



6 



A. DEMLOW, J. GUZMAN, AND A. SCHATZ 



satisfies IMI^/j f Gl \ < Cdr d ' , < j < r. We may then apply the assumptions Al 
through A3 to the pairs G\ and G2, and G2 and G. 
Using (|3.3p . we first compute that 

(3.10) \\ u h\\ 2 H i(G ) < WuuhWnuG) < CL(uju h ,u)u h ). 

Using the fact that HVwHjwq) < §, we compute that for any e > 0, 

L(u;u h ,u)Uh) = L(u hl uj 2 u h ) 

— / Uh[AV(ujUh)Vu + UhA\7uj\7uj + AVu;V(ujUh) + u>Uhb\7u)] dx 



(3.11) 



<\L(u h ,uj 2 u h )\ + C-^\\u h \\l 2(G) +e||ww/,|||- l(G) . 



Next we use (JHHJ), ([3~7| , and the fact that \\L0 2 u h \\m(G) < ||^Uft ||ffi(G) + §ll u ft||L 2 (G) 
to compute 

L(u h , Lo 2 u h ) =L(u h , tu 2 u h - I(u) 2 u h )) 

^ 3 ' 12 ^ <C h T\\Uh\\H^T)(j\uU h \ H i m + -p\\u h \\L 2 (T))- 

Using (|3.6|) and the fact that < 1, we have for e as above that 

Ch T \\u h \\ H i {T) (huju h \ H i {T) + -L||u ft ||i, 2 (T)) 

(3.13) C 

<-^p\\uh\\ 2 L2{ T) +t\uJu h \ 2 H1{T) . 

Inserting ([3~T3]) into (f3~T2]) . noting that T n G 2 ^ implies that T C G (since 
max TnG7 i0 hj< < j ) and carrying out further elementary manipulations then yields 
that for e > 0, 

C 

(3.14) L{u hl uj 2 u h ) < —p |K||! 2(G) + e||o;ufc||^i (G) . 

Inserting (13.141 into p. lip and the result into (|3.10|) yields 

G 

(3-15) llwfclll-i(G) < ^2 H^llL(G) + 2e||wu h ||^i (G) . 

Taking e = ^ so that we may kick back the last term above, employing the triangle 
inequality, and inserting the result into (|3.10[) then completes the proof of (|3.9[) . □ 

We now prove a local energy error estimate. In our proof below we shall follow 
NS74J by using a local finite element projection in order to split the finite element 
error into an approximation error and a "discrete harmonic" term which may be 
bounded using Lemma [3~3l We note, however, that the use of a local finite element 
projection is not necessary, and our final local error estimate may in fact be proved 
with some simple modifications to the proof of Lemma |3~31 above. These two styles 
of proof are essentially equivalent. Local finite element projections have been used 
for example in |NS74| , |SW77j , [SW95] , and |AL95| in order to prove local a priori 
error estimates. The methodology of Lemma 13.31 in which no local projections are 
used has been employed in for example [Dem04 and Guz06j in order to prove 
local a priori error estimates and in LN03 and [Dem07 in order to prove local a 
posteriori error estimates. 
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Theorem 3.4. Let Go C G C ilo be given, and let dist(Go,dG \ dil) — d with 
d < min{2do, di} where do is the parameter defined in the assumption Rl and d\ 
is defined in Remark \3.1[ Let also Al, A2, and A3 hold with d = j^. In addition 
let max TnG ^0 < jg. Then 

\\u-u h \\ H i {Go) <C min (\\u-x\\hhg) + k\\u ~ x\\ l 2 (g)) 
(3.16) i 

+ C-\\u-u h \\ L2{G) . 



Here C depends only on the constant C in S2.2\) and the coefficients of L. 

Proof. We assume that Go is the intersection of a ball Bd of radius | with fl ; the 
general case may be proved using a covering argument. Let Gi be the intersection 
with Qo of a ball having the same center as Go and having radius ^ , and without 
loss of generality let G be the corresponding ball of radius d. Let then ui G Gq 30 (G) 
be a cutoff function which is 1 on G\ and which satisfies IMIw^ ( G ) ^ Cdr^ , 
< J < r. Note that we may apply Lemma \3. 31 with Go on the left hand side of 
the estimate (|3 .9[) and Gi on the right hand side. 

Next we let P{u>u) be a local finite element projection of uju. In particular, we 
let P(lju) £ S D ("1 #<(G) with u h - P(lju) = on T D n dG ± satisfy 

(3.17) L{um - ?H, v h ) = 0, v h e S D ,o n H< (G). 

The local coercivity condition f|3. 3|) then implies the stability estimate 



(3.18) ||PM|[ K i (G) <C\\uju\\ h i (g) . 

Recalling that Uh — P(wu) = on Tr> H <9Gi while employing (|3.9[) and using 
(|3.4p while recalling that oj = 1 on Gi, we compute that 

||W - U/i||ffi( Go ) < - - p ( ww )llff 1 (G ) + ll p ( ww ) ~ ^IIhHGo) 

G 

<||wu - P(ww)|| H i (G) + —\\P(uju) - Ufe|| i2 ( Gl ) 

(3.19) c* 

<||wu - P(ww)||iri(G) + --j(||P(wu) - uju\\ L2[Gi) + |]u - u h \\ L2(Gl) ) 

C 

<C\\uju - P(uiu)\\ H i (G) + u h \\ L2{Gl) . 

Next we employing the triangle inequality along with (|3 . 1 8|) while recalling that 
IMIw^ (G 2 ) — Cd~ J in order to find that 

\\uju- P(uju)\\ h i {g) <C\\uju\\ h i {g) 

(3.20) 1 

<C(\\u\\ HHG) + -\\u\\ L2iG) ). 

In order to complete the proof of (|3.16p , we first insert (|3.20p into (|3.19|) and 
finally write u - uh = (u — x) + (x ~ u h) with S D ,o- 

□ 
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